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ABSTRACT 

The Galform semi-analytic model of galaxy formation is used to explore the 
mechanisms primarily responsible for the three types of galaxies seen in the local 
universe: bulge, bulge+disk and disk, identified with the visual morphological types 
E, SO/a-Sbc, and Sc-Scd, respectively. With a suitable choice of parameters the Gal- 
form model can accurately reproduce the observed local if s -band luminosity function 
(LF) for galaxies split by visual morphological type. The successful set of model pa- 
rameters is used to populate the Millennium Simulation with 9.4 million galaxies and 
their dark matter halos. The resulting catalogue is then used to explore the evolution 
of galaxies through cosmic history. The model predictions concur with recent obser- 
vational results including the galaxy merger rate, the star formation rate and the 
seemingly anti-hierarchical evolution of ellipticals. However, the model also predicts 
significant evolution of the elliptical galaxy LF that is not observed. The discrepancy 
raises the possibility that samples of z <~ 1 galaxies which have been selected using 
colour and morphological criteria may be contaminated with galaxies that are not 
actually ellipticals. 

Key words: galaxies: elliptical and lenticular, cD; galaxies: evolution; galaxies: for- 
mation; galaxies: luminosity function; galaxies: spiral; galaxies: structure 
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1 INTRODUCTION 

The current paradigm, constrained in part by the K, 
(hereafter A'-baiid) luminosity function (LF) (e.g 



-band 



Benson 



|et al. 2003 1 , has galaxy disks forming through a combina- 
tion o f cold gas accretion (|Weinberg et al.||2004| ) and feed- 
back (Oppenheimer & Dave 20061, with bulges resulting 



from mergers (e.g., Barnes & Hernquist 1992; Hopkins et al 
|2008| |Masjedi et al.||2008[ ). Thus, the diversity of galaxy 
morphologies seen today represents the culmination of mul- 
tiple evolutionary paths. These end points in galaxy evolu- 
tion are captured in a taxonomy devised b y|Hubble| (|1936[) 



and refined by de Vaucouleurs (19591 and de Vaucouleurs 



et al. ( 1991 1 that is based on the relative prominence of 



the stellar bulge and the degree of resolution of the spiral 
arms. Despite attempts to emulate this classification scheme 



using quantitative proxies such as colours (e.g., Strateva 
et al.|200Tj ), concentration-asymmetry in dices (e.g., |Watan- 



abe et al. 1985 Abraham et al. 1996 Bell et al. 20031, 



the Gini-M20 classifiers (Lotz et al. 20041 and parametric 
methods involving bulge/disk decompositions (jSchade et al.| 



1996| |Ratnatunga et aT|l 999 Sim ard et al.|2002[ ) there are 
no quantitative measures that uniquely segregate galaxies 
into the morphological groups represented by the de Vau- 
couleurs numerical T system (e.g., Graham fc Worley|2 008 ) . 



Thus, proxies are merely a useful stopgap until the difficult 
and time consuming task of assigning visible morphologies 
to the ever-growing number of cataloged galaxies can be ac- 



complisheCj 

Fortunately, the time-consuming task of assigning visi- 
ble morphologies has now been completed for the vast ma- 
jority of nearby galaxies. The principal aim of this paper 



is therefore to use the K-band LF for nearby galaxies (De- 



vereux et al. 2009 1 to constrain models of their formation 



and evolution. The Tf-band LF is of particular interest be- 
cause of its relevance to understanding galaxy evolution in 
the context of Lambda Cold Dark Matter (ACDM) cosmol- 
ogy; at z — 0, the if-band LF traces the stellar mass accu- 
mulated in galaxies at a wavelength where interstellar ex- 



tinction is minimal (Devereux et al. 1987 Bell & de Jong 
[2000} |Bell et al.||2003[ ). Additionally, the functional forms 
for the if-band LFs distinguish between bulge dominated 
and disk dominated systems which suggests that at least 
two quite distinct galaxy formation mechanisms are at work 
to produce the diversity of morphological types seen in the 
local universe. Previous attempts to model the LFs in the 



context of hierarchical clustering scenarios (e.g., Cole et al 



1 See the Galaxy Zoo Project, http://www.galaxyzoo.org/ 
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Morphological 




B/T distribution 


class 


B/T Range 


20% 


50% 80% 


E 


0.92 < _B/T K sS 1.00 


0.99 


1.00 1.00 


SO-Sbc 


0.11 < B/T K < 0.92 


0.22 


0.39 0.65 


Sc-Scd 


0.00 < B/T K < 0.11 


0.00 


0.01 0.05 



2 MODEL 

The Galform semi-analytic model has been described most 



Table 1. The Table presents the mapping between morphologi- 
cal type and K-band bulge-to-total ratio required for our best-fit 
model to reproduce the relative abundances of each class in the 
interval —23.5 < Mk — 51og 10 /i —23. Additionally, the final 
three columns show the percentage of galaxies in each class with 
Mk — 5 log h < — 21 and B/T ratios less than the value indicated. 
For example, 80% of bright SO-Sbc galaxies have a B/T less than 
0.65. 



|2000| |Benson et al.|2003| have focussed on the Universal or 
Total iif-band LF, which is the sum over all galaxy types, 
because they are based on samples for which there are no 
pre-existing morphological assignments. In this paper, the 
_K"-band LFs are modeled for galaxies segregated by visible 
morphology in order to better understand: 1) the origin and 
evolution of bulgeless disk galaxies, represented in the ob- 
servable universe by Sc-Scd galaxies; 2) the role of secular 
processes in building the bulges of spiral (SO/a — Sbc) and 
lenticular galaxies and; 3) the galaxy merger history that 
led to the formation of elliptical galaxies. Our analysis dif- 



fers from |Parry et al.| ( |2009| > in that the GALFORM model 
predictions are constrained by the observed A"-band LFs for 
galaxies of different morphological type. Those model pre- 
dictions are then compared with recent observational results 
including the type averaged galaxy merger rate and the type 
specific star formation rate, both quantified as a function of 
redshift. 

Observationally, there is a weak dependence of K-band 
B /T on visual morphology ( |Graham fc Worlcy 2008 ) which 
is used in this paper to segregate the model galaxies into 
three broad groups representing the three types of galax- 
ies seen in the local universe, namely bulge, bulge+disk and 
diskQ These three groups are identified with the visual mor- 
phological types E, SO/a-Sbc, and Sc-Scd, respectively. Fol- 
lowing [GrahanT^&^Vbr ley ( 2008]) , these three broad groups 
are identified in the models by adopting the range of B/T 
ratio listed in Table [T] The table also indicates the distri- 
bution of B/T within each morphological class, illustrating 
that most disk galaxies (including SO's) have a rather low 
B/T ratio. For example, 80% of SO and later type galaxies 
brighter than M K - 5 log 10 h = -21 have B/T ^ 0.54. If the 
magnitude limit is lowered to Mk — 5 log 10 h = —19, we find 
that 80% of such galaxies have B/T ^ 0.35. This is in agree- 
ment with the broad conclusion of |Graham fc Worley| ( |2008[ ) 
that most disk galaxies (including SO's) have low B/T ratios 
(e.g. B/T <l/3). 



2 In principle, other outputs from the model could be used to 
provide additional information on morphological type. However 



|Graham &: Worlcy (2008) show that bulge/disk flux ratio is the 
most useful discriminant as it varies more, as a function of mor- 
phological type, than other indicators, such as bulge/disk size 
ratio, for example. 



recently by Bower et al. ( 2006 1 and the reader is referred 
to that paper for a full description. In the model, galax- 
ies consist of an exponential disk containing both gas and 
stars and a spheroidal stellar component — the morphology 
of the model galaxies is therefore determined by the rela- 
tive amount of light in each component. A few key points 
regarding morphological evolution of model galaxies bear 
reiterating however: 

• Galaxies initially form through the cooling of gas from 
their surrounding dark matter halo. This gas is assumed 
to conserve its angular momentum during collapse and to 
therefore form a disk. As such, a galaxy formed in isolation 
will always be a pure disk galaxy unless it is sufficiently 
self-gravitating to become unstable. 

• Disks which become too self-gravitating become unsta- 
ble to the formation of global perturbations and therefore 
destroy themselves. Instability is judged using the criterion 



proposed by Efstathiou et al. ( 1982 1. For a disk to be stable 
it must satisfy 



I'D 



(GMv/Rv) 1 /* > (1) 

where v is the circular velocity of the disk at its half-mass 
radius, Ru and Mo is the mass of the disk. The value of 
e for purely stellar disks was found to be 1.1 by Efstathiou| 
et al. ( 1982 I, although lower values may be appropriate for 
gas rich disks ( Christodoulou et al.|[T995 |. We treat e as a 
free parameter. If a disk is deemed to be unstable according 
to this criterion, the disk is destroyed, with any stars and 
gas present forming a spheroidal system in which the gas is 
then rapidly turned into stari|^] The end result is a purely 
spheroidal (a.k.a. elliptical) system. 

• Mergers between galaxies can also create a spheroidal 
component. Galaxy-galaxy mergers are driven by dynamical 
friction on a satellite galaxy (and its halo) orbiting in a 
larger dark matter halo. Merging timescales are computed 
using a modified form of Chandrasekhar's dynamical friction 
equation (Lacey & Cole|1993|. The timescales are scaled by 



.,,„,, a parameter of the model, to allow some adjustment 
of merger rates. Mergers are split into two categories: 

- Minor Mergers: These are mergers with M2/M1 < 
/eiiip where Mi and M2 < Mi are the baryonic masses 
of the galaxies involved in the merger and / c in p is a pa- 
rameter of the model (see below). In a minor merger, the 
stars from the secondary (lower mass) galaxy are added 
to the spheroidal component of the primary (more mas- 
sive) galaxy while any gas from the secondary is added 
to the disk of the primary. A galaxy can therefore grow 
a spheroidal component through the cumulative effects of 
many minor mergers. 

- Major Mergers: These are mergers with M2/M1 ^ 
/eiiip and cause the destruction of any prexisting stellar 
disks, such that all stars from the merging galaxies become 



3 Since the treatment of this instability process in semi-analytic 
models is somewhat uncertain, an alternative and less dramatic 
implementation is presented in the Appendix |A2| The results do 
not change the qualitative conclusions of this paper. 
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part of the spheroidal component in the merger remnant. 
The fate of any gas in the merging galaxies also depends 
on the merger mass ratio: 

- If M2/M1 < /burst, where /burst is a parameter of the 
model, then the merger does not trigger a burst of star 
formation and any gas forms a disk in the merger remnant 
(and may subsequently form stars in that disk). 

- If M2 /Mi ^ /burst then the merger triggers a burst of 
star formation and any gas present in the merging galaxies 
is turned into stars on a short timescale and becomes part 
of the spheroidal component of the merger remnant. 

• Once a spheroidal component has formed, subsequent 
gas cooling can cause a disk to regrow around the spheroid, 
resulting in a galaxy with intermediate bulge-to-total ratio. 

• The model includes supernovae feedback which is as- 
sumed to drive gas out of a galaxy at a rate propor- 
tional to the star formation rate: M out = /3M+, where 
P = (Vhot /V) ahot , V is the circular speed of the galaxy 
and Vhot and « no t are parameters of the model. Vhot may 
differ for quiescent star formation in disks and bursts of star 
formation. 



The |Bower et~aLl |2006 ! ) model was not constructed us- 
ing morphological data as a constraint. Our initial goal, 
therefore, is to adjust the parameters of the |Bower et al.| 
(2006) model so that it reproduces the observed morpho- 



logically segregated LFs of Devereux et al. (2009), shown in 



Figure [T] A search for an alternative set of values for the six 
model parameters that match the observed morphologically 
selected LFs yielded the values listed in Table [2] The best 
fit model LFs (lines) are compared with the observed LFs 
from |Devereux et al.| ( |2009[ ) (points) in Figure [I] Morpho- 
logical class is indicated by colour in the figure. The relative 
abundance of each type in the model is forced to agree with 
the data in the interval —23.5 < Mk — 51og 10 /i ^ —23.0 
(this is how the mapping between B/T and morphological 
class listed in Table [I] was determined) but the morpholog- 
ical mix is not enforced outside this magnitude range. Also 
shown, for reference, are the results from the unmodified 



Bower et al. (2006) model (dashed lines). The comparison 
shows that the model results are quite similar. 

The model LFs for ellipticals, lenticulars and bulge- 
dominated spirals are peaked and decline toward both higher 
and lower luminosities. In contrast, the Sc-Scd class contin- 
ues to rise towards low luminosities and dominates the total 
LF for A/k — 5 log 10 h > —22. The black line and points indi- 
cate the total (i.e. all morphological types) LF. Not surpris- 
ingly, the model prediction here is also in good agreement 
with the data. This is to be expected however, since the 
model parameters are optimised to yield a satisfactory rep- 
resentation of the observed LFs for each of the three broad 
morphological types and the total LF is just the sum of 



these. Additionally, Bower et al. ( 2006 1 have already demon- 



strated that the model can fit the observed total K-band 
LF. The parameters yielding the best-fit model (Table |5| 



are similar to those adopted by Bower et al. (20061: disk 
instability estimates are unchanged, but the rate of galaxy 
merging has been reduced by a factor of 2, a lower mass ratio 
in mergers is required to disrupt stellar disks, a higher ratio 
is required to trigger a burst, while the feedback in bursts, 
but not disks, has been weakened. 




Figure 1. Morphologically selected K-band LFs for the best fit 
model (Table 2) are shown by solid lines and are compared to 
the observed LFs of |Devereux et ah] ( |2009| ) shown by the points. 
Colours indicate morphological class as indicated by the key in 
the figure. Galaxies are split into three morphological classes ac- 
cording to Table 1. For reference, dashed lines show results from 
the |Bower et al!] l |2006[ l model (i.e. without the parameter modi- 
fications introduced in this work) using the same cuts on B/T to 
define morphological classes. 



Table 2. Parameters of the best-fit model. 



Value 



Parameter Best Fit Bower et al. 2006 



c 


0.8 


0.8 


TOmrg 


3.0 


1.5 


/ellip 


0.15 


0.3 


/burst 


0.2 


0.1 


^hot,disk 


485km/s 


485km/s 


^hot, burst 


300km/s 


485km/s 



3 RESULTS 

Having found a successful set of model parameters, they are 
used to populate the Millennium Simulation ( Springcl ct al. 
2005) with galaxies. To do this, dark matter halo merger 



trees are extracted from the Millennium Simulation using 
the techniques described by [Harker et al.| ([2006 ) . The trees 
are then populated with galaxies in the manner described by 



Bower et al. (20061, but using the model parameters found 
in the previous section. This results in a catalogue of 9.4 
million galaxies and their dark matter halos. The catalogue 
provides the magnitudes, morphologies and other physical 
characteristics for all of the galaxies, including their evolu- 
tion through cosmic history. The catalogue can therefore be 
used to explore galaxy evolution. Our goal is to reveal the 
key physical processes at work that produce the three broad 
types of galaxies in the observable universe, namely bulges, 
bulges+disks and disks. The resolution of the Millennium 
Simulation, as in any N-body simulation, is limited. In Ap- 
pendix [Alia study of the effects of resolution demonstrates 
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M K -5log l0 h M K -5log, h 



Figure 2. The median halo mass in which galaxies formed (solid 
lines) and the mass of their current halo at z = (dashed lines) as 
a function of K-band luminosity and morphological class (colour 
coding as in Fig. The current halo mass is the same as the 
mass of the halo in which the galaxy formed for central galaxies 
but may be significantly larger for satellite galaxies. 



that the morphological features of the model galaxies are 
well converged at the resolution of the Millennium Simula- 
tion. 



3.1 Trends with Morphological Class 

3.1.1 Dark Matter Halo Properties 

The relationship between galaxy luminosity and dark matter 
halo mass is illustrated in Fig. [2] as a function of morpho- 
logical type. The solid lines in Fig. [2] indicate the mass of 
the halo in which a galaxy formed, defined here following 



Cole et al. ( 2000[ ) such that halos are labelled as "reform- 
ing" each time they undergo a doubling of their mass. As 
may be expected there is a correlation between galaxy lumi- 
nosity and the median mass of the dark matter halo in which 
such galaxies formed. The slope of the relation is quite flat 
however, showing a range of a factor of ~20 in halo mass 
over a factor of ~100 range in luminosity. 

There is little dependence of median formation halo 
mass on morphological type, except for the E class which 
form in halos around a factor of ~ 3 more massive for a 
given luminosity. This may be expected for the E class if such 
galaxies are built up via merging rather than continuous gas 
accretion. To attain the same luminosity as its accretion- 
fueled counterpart, a merger-built galaxy must form in a 
more massive halo so that the rate of merging will be en- 
hanced. 

The dashed lines in Fig. [2] indicate the mass of the dark 
matter halo in which a galaxy lives at z — 0. In the case of 
satellite galaxies, this is the mass of the isolated dark matter 
halo in which they orbit, not the mass of the subhalo within 
which they reside. The difference between the dashed and 
solid lines is therefore indicative of how much mass has accu- 
mulated in the dark matter halo since galaxies formed. Dark 



Figure 3. The Sc-Scd K-band LF (blue line and points) versus 
the median formation halo mass function (black line) scaled by a 
M/L ratio. 



matter halos may grow through mergers of satellites and by 
accretion. With the exception of Sc-Scd galaxies, the mass 
of the halo in which galaxies live at z — greatly exceeds 
the mass of halo in which they formed, particularly for the 
lower luminosity systems (a simple consequence of the fact 
that the low mass halos in which these faint galaxies form 
are highly likely to be subsumed into much more massive 
halos by z = by becoming satellites of those halos). 

Interestingly, Sc-Scd galaxies of all luminosities occupy 
dark matter halos that are only about a factor of two more 
massive than those in which they formed. Thus, disk dom- 
inated galaxies must have avoided major mergers by virtue 
of living in dark matter halos having rather quiet merging 
histories. There are a wide variety of possible merger histo- 
ries for halos of a given mass ( |Stewart et al.| 2008) including 
some fraction which will have avoided major mergers. It is 
known that halos in low density environments experience 
lower merger rates than their counterparts in high density 
environments ( Fakhouri fc Ma|2009 1 and so are more likely 
to have a quiet merger history. Thus, Sc-Scd type galaxies 
are distinguished from the other morphological types hav- 
ing formed preferentially in low density environments. 65% 
of Sc-Scd galaxies brighter than Mr — 5 log 10 h — —20 are 
central galaxies, ie. they are the most massive galaxy in their 
halo, and they grow primarily through gas accretion and not 
mergers. Consequently, the Sc-Scd K-band luminosity func- 
tion has a similar slope to the median formation halo mass 
function as illustrated in Fig. [3] 



3.1.2 Tvmescales 



While interesting from a theoretical standpoint, dark mat- 
ter halo masses cannot be easily measured for most galax- 
ies. (Although weak lensing ( Mandclb aum et al.||2006 l and 
clustering ( Mandelbaum et al.|2009[ ) measures can constrain 
halo mass for classes of galaxies.) Greater utility is provided 




Figure 4. The stellar mass-weighted age of galaxies as a function 
of K-band magnitude and morphological class (colour coding as 
in Fig.Q. Solid lines indicate the median ages for galaxies in each 
morphological class while error bars indicate the ±1<t interval of 
the distribution of ages. 



by examining observable quantities more directly related to 
the galaxies themselves, beginning with timescales. 

Figure [4] shows the median stellar mass- weighted age of 
galaxies, defined as 



(tage) 



(2) 



where to is the current age of the Universe and p*(t) is the 
star formation rate in a galaxy at time t, split by morpho- 
logical type as a function of their K-band magnitude. Disk 
dominated galaxies (blue line) are distinguished once again, 
this time by hosting a uniformly much younger stellar pop- 
ulation than any other morphological type. This suggests 
that disk galaxies of all luminosities have been continuously 
forming stars since their formation which the model reveals 
is at a rate that is proportional to the baryonic mass of the 
galaxy. In contrast the ellipticals, lenticulars and bulge dom- 
inated spirals contain much older stellar populations with 
median ages of lOGyr. This does not mean that the galaxies 
formed lOGyr ago, merely that the stars now dominating the 
baryonic mass of the galaxies typically formed at that time. 
As we shall see in |3.1.4| massive ellipticals were actually 
assembled much more recently. 



3.1.3 Bursts and Bulge Formation 

Figure [5] examines the fraction of a galaxy's bulge mass that 
was made in the most recent burst of star formation (solid 
lines) and the fraction formed through disk instabilities (the 
remainder being made through minor or major mergers). 
Interestingly, ~30-50% of the bulge mass in lenticulars and 
bulge dominated spirals (SO/a - Sbc) resulted from disk in- 
stabilities with the remainder created in the most recent 
burst of star formation, which occurrred ^lOGyr ago, that 
itself could have been triggered by a disk instability or a 



Figure 5. The fraction of a galaxy's bulge mass which was made 
in the most recent burst of star formation (solid lines) and disk 
instabilities (dashed lines) as a function of K-band magnitude for 
each morphological class (colour coding as in Fig. [Tjl. 



merger. Unfortunately, there is no way to distinguish be- 
tween the two possibilties in this figure. The importance of 
these processes is greatly diminished in elliptical galaxies, 
because they do not possess diskfQ and in disk dominated 
Sc-Scd galaxies, because they have small or non existent 
bulgefj^] 

3.1.4 Formation and Assembly Redshift Distributions 

Figure [6] shows the distribution of star formation redshifts, 
z*, and assembly redshifts, z asscm . The star formation red- 
shift for a galaxy is defined to be the redshift at which half 
of the final stellar mass of the galaxy has been formed. The 
assembly redshift is the redshift at which half of the final 
stellar mass has been assembled into a single progenitor. 
Clearly for a given galaxy 2 assC m < z*- The most striking 



4 The dashed lines reflect the mass of stars that were formed 
through disk instabilities in all progenitors. Thus the dashed red- 
line indicates that disk instabilties may have been an important 
mechanism for the progenitors of todays ellipticals. We note that 
the fraction of an elliptical's mass which formed as a result of disk 
instability events is sensitive to the details of the treatment of 
such instabilities in our model. With our standard treatment, on 
averge 29% of the spheroid mass in bright (Mk — 5 log 1( ) h < —23) 
ellipticals is made via disk instability events. Using a more mod- 
erate treatment of instability events in which a minimal amount 
of disk is converted into spheroid (see Appendix |A2| this num- 
ber is reduced to 9%. These fractions are relatively insensitive to 
the resolution of the merger trees. The resolution study described 
in Appendix |A1| indicates that they change by only 3% for the 
brightest ellipticals as the resolution of the merger trees is in- 
creased by a factor of 16. For the faintest ellipticals shown (which 
are closer to the resolution limit of our calculations) the fraction 
changes by only 7% with the same increase in resolution. 

5 However, the models do indicate that disk instabilities may 
be an important bulge formation mechanism in the luminous, 
Mk — 51og 10 h < — 23 mag, Sc-Scd galaxies. 
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2 4 2 4 6 

z 

Figure 6. The normalized distributions of (the redshift at which half of a galaxy's final number of stars were formed; dashed lines) 
and z a ssem (the redshift at which half of the galaxy's final stellar mass was assembled into a single progenitor; solid lines) are shown split 
by morphological class (colour coding as in Fig.JTJ. Each panel corresponds to a different K-band magnitude interval as indicated in the 
panel. The distribution if normalized such that the integral under the curve is unity. 



feature of these distributions is that the most luminous -25 
^ Mr — 5 log 10 h ^ -24 elliptical galaxies have a median 
star formation redshift z* ~ 3 while the median assembly 
redshift is much lower z assom sC 3. This clearly shows that 
such galaxies are formed almost entirely through the merg- 
ing of pre-existing systems that formed at higher redshift. 
The difference between star formation and assembly times 
diminishes for lower luminosity ellipticals, suggesting fewer 
mergers. Also, the stellar population is younger in lower lu- 
minosity ellipticals compared to their more luminous coun- 
terparts by ~ 1 Gyr. Since the other morphological types 
(SO-Sbc) show remarkably little difference between their star 
formation and assembly redshift distributions, it is primarily 
the luminous ellipticals, -25 ^ Mk — 51og 10 h ^ -24 whose 
growth is dominated by the delivery of pre-formed stars via 
mergers that happened recently, since z ~ 3. 

The star formation and assembly redshifts are virtually 
identical for lenticular and spiral galaxies (SO-Scd) which 
suggests that major mergers do not play an important role 
in the evolution of these galaxy types. The z* distributions 
are strongly peaked at low redshifts for the disk dominated 
Sc-Scd galaxies. Thus, the models indicate that they should 
contain the youngest stellar populations with ages ^ 6 Gyrs. 
On the other hand, the z* distributions for lenticular and 
bulge dominated spirals become broader with decreasing lu- 



minosity suggesting a heterogeneous group with stellar pop- 
ulations spanning a wide range in age. 



4 DISCUSSION 

In order to understand both the trends found in ^3] and 
the morphological evolutionary paths of galaxies, we exam- 
ine formation tree diagrams such as Fig. [7] These show the 
detailed formation history of individual galaxies and there- 
fore provide a more detailed examination of morphological 
evolution. The tree shown in Fig. [7] is for a luminous ellip- 
tical galaxy. As expected from the results of ^3] this galaxy 
shows significant merging activity. There is a clearly identi- 
fiable main progenitor which has been a pure elliptical since 
z « 0.5. Half of the stars in the final galaxy were formed 
by z = 3.75 but where not assembled into a single galaxy 
until z = 0.26. Approximately one quarter of the stars in 
the galaxy were formed as a result of a burst triggered by a 
disk instability event which occurred at z ~ 5. 

In Figs. H] through [10] we show tree diagrams for ran- 
domly selected galaxies in the E, SO/a-Sbc and Sc-Scd 
classes for a range of luminosities. (Note that the vertical 
axis of each plot is adjusted to stretch from z = to the 
earliest redshift at which the galaxy has an identifiable pro- 
genitor.) The E class shows a clear transition from forming 
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M,(-51og(h) = -25.5 
B/T K =1.00 (EJ 




0.5 



M/M„ 



Figure 7. A formation tree diagram for a single galaxy identified 
at z = 0. The K-band magnitude and bulge-to-total ratio of the 
galaxy at z = are shown in the upper right-hand corner. The 
y-axis plots redshift, while the x-axis plots stellar mass relative to 
that of the 2 = galaxy. The small pie-charts indicate the final 
galaxy and any progenitor galaxies (shown down to progenitors 
with a mass of 0.1% of the final galaxy). Black lines connecting 
galaxies indicate mergers between galaxies. The radius of each 
pie-chart is proportional to the cube root of the stellar mass of 
the galaxy, while the segments of the pie chart show the relative 
fractions of stars in the galaxy's disk (blue), and bulge (red for 
stars formed by or brought in by mergers, green for stars formed 
or brought in by disk instabilities). The dotted circles around each 
galaxy are proportional to (Q\ > /Qo)My l and so will coincide with 
the size of the galaxy pie chart if a galaxy has succeeded in con- 
verting all available baryons in its halo into stars. Coloured lines 
indicate the total stellar mass (magenta) as well as that in the 
disk (blue), bulge [mergers] (red) and bulge [instabilities] (green) 
as a function of redshift. Solid lines show masses in the most 
massive progenitor, while dashed lines show the mass summed 
over all progenitors. 



via multiple mergers in the high-luminosity regime to form- 
ing via disk instabilities in the low luminosity regime. (A 
disk instabiity event in these diagrams is apparent as the 
solid [main progenitor] blue [disk] line drops to zero while the 
solid green and red lines jump to larger values.) These low 
luminosity ellipticals start out as disks, gradually converting 
their gas into stars. Eventually, the disk becomes gravita- 
tionally unstable, triggering a burst and converting all stars 
and gas into a spheroid. The galaxy evolves passively after 
this time. 

It is apparent that the number of progenitors is greatly 



reduced for the lower luminosity galaxies^] This is apparent 
for all morphological types except for the Sc-Scd class and 
is a consequence of the shape of the galaxy mass function in 
systems of differing halo mass (e.g. Benson et al.|2003 l — for 
the most massive halos the galaxy mass function is a very 
steep function of halo mass, such that there are many more 
galaxies of slightly lower mass with which the central galaxy 
can merge. For lower mass halos, the galaxy mass function 
is shallower, yielding relatively fewer galaxies of lower mass. 
This results in fewer galaxies with which to merge and fewer 
progenitors. The shape of the galaxy mass function in our 
model arises because of supernovae feedback which is rela- 
tively stronger in lower mass halos. 

The SO/a-Sbc class has a much quieter history. Only the 
most luminous examples Mk — 5 log 10 — 25 are involved 
in a few significant mergers. Lower luminosity systems expe- 
rience periods of early growth characterized by multiple disk 
instability events that cause the disk mass to rise and then 
fall sharply. More recent growth, since z~2, occurs when 
the disk has stabilized and is able to form stars quiescently 
until the gas supply is exhausted, typically between z = 1 
and z — 2 for the lower luminosity systems and z~0 for the 
luminous systems. 

The Sc-Scd class has the most uneventful formation his- 
tory of all. This is, of course, obvious from the fact that 
this class is defined to have almost no spheroidal compo- 
nent. These galaxies must therefore have experienced no 
significant instability or merger events. The tree diagrams 
of Fig. [To] show that this type of galaxy is characterized by 
steady star formation in a disk fueled by cosmic infall, which 
continues either until z = 0, or until the galaxy becomes a 
satellite and therefore loses its gas supply. 

4.1 The Merger History as a Function of 
Morphological Type 

Figure [TT] illustrates the merger history as a function of mor- 
phological type. F(z) is the fraction of galaxies with mass 
>2.5 x 10 10 M Q (corresponding to M K - 51og 10 - 22.7) 
that have undergone a major merger (M2/Mi>0.25) within 
the past lGyr. These selection criteria allow a direct compar- 
ison with the observationally determined rates of |Jogee et al.| 
(2009) and |L6pez-Sanjuan et alT| ( |2009[ ) . Galaxies are segre- 
gated according to their z = morphology and the reported 
merger fractions resulted from integrating their respective 
merger trees over z. The model predicts an average merger 
fraction integrated over all galaxy types of ~ 2% which com- 
pares favorably with the observed lower limit to the major 
merger fraction for the sample SI of |Jogee et al.| ( [2009[ ) and 
the merger fraction measured at z = 0.6 for a slightly more 



luminous sample by Lopez-Sanjuan et al. (2009 1. The model 
reveals that the merger rate is highest for elliptical galaxies; 
~ 8% of ellipticals at z — resulted from a major merger 

e Galaxies of luminosity similar to that in panel d of this figure 
are affected by the resolution limit of the simulation as the halos 
that they occupy typically contain around 50 particles, preventing 
their lower mass progenitors from being resolved. For galaxies of 
luminosity comparable to that in panel c (and brighter galaxies) 
resolution is not an issue — they inhabit halos typically containing 
500 particles, sufficient to permit many lower mass progenitor 
halos to be resolved in the simulation. 
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that occurred in the past 1 Gyr. The model elliptical merger 
rate is significantly higher than the merger rate estimated by 



Masjedi et al. (2008) for a sample of luminous red galaxies 
in close pairs. The model elliptical merger rate predicts sig- 
nificant evolution of the elliptical LF, as discussed in more 
detail in the next section. Very few of the z — lenticular 
and (SO/a-Sbc) spiral galaxies resulted from a major merger 
and none of the Sc-Scd disk galaxies. Thus, the model pre- 
dicts that Sc-Scd galaxies exist in the local universe because 
they have not been involved in any collisions. 

We compare our predicted merger rates with other re- 
cent attempts to measure this quantity from either N-body 
simulations of dark matter or from semi-analytic models of 
galaxy formation. Before doing so, we remind the reader that 
our merger rate is that of galaxies, and so there is no reason 
to expect it to coincide with rates of, for example, dark mat- 
ter halo merging or dark matter subhalo destruction (both 
of which are commonly measured quantities). In particular, 



corrections of approximately k = —0.75 and e = 0.55, to be 
compared to k — —1.05 ± 0.15 and e = 0.4 from Bel l et al.| 



Stewart et al. (20091 examined the rates of both halo merg- 



ing (when a halo first crosses the virial radius of a larger 
halo, thereby becoming a subhalo) and subhalo destruction 
(when a subhalo loses 90% of its mass) in an N-body sim- 
ulation of dark matter. Stewart et al. ( 2009 1 find merger 



rates which increase monotonically with increasing redshift, 
in disagreement with our own findings. However, we caution 
that these results are based purely on dark matter properties 
(with a simple, redshift dependent mapping of galaxy lumi- 
nosity to dark matter halo mass) and so do not measure the 
same quantity as presented in this work. Similar results were 
found by Genel et al.| ( |2009| ). On the other hand, Bertone & 
Conselice ( 2009 1 find merger fractions comparable to ours, 
with a broad peak around z = 1.5 and a gradual decline 
towards higher redshifts. The agreement is not surprising, 
however, as they derive merger rates from a semi-analytic 
model similar to the one used here. 



4.2 The Evolution of the Elliptical Galaxy 
Luminosity Function 

The model predicts significant evolution of the elliptical 
galaxy K-band luminosity function as illustrated in Fig- 
ure !12l The LF at each redshift has been k-corrected and fur- 
ther corrected for passive evolution to z — so that the LFs 
may be compared on an equal basis. Within the Galform 
model, k-corrections and passive evolution can be computed 
precisely. The model predicts the full SED of each galaxy 
and so we can simply shift the filter to the galaxy rest frame 
and compute the absolute magnitude which will then au- 
tomatically include the k-correction. For passive evolution, 
the model predicts the mix of stars in each galaxy at each 
redshift, including their distribution of ages and metallic- 
ities. To compute the magnitude of the galaxy at z — 
assuming passive evolution (i.e. no further star formation) 
we simply artificially increase the ages of these stars by the 
lookback time to the current redshift and then sum their 
SEDs to give the net SED of the galaxy passively evolved 
to z = 0. The passively evolved magnitude of the galaxy 
is then trivially found by integrating under the appropri- 
ate filter response function. Our k and evolution corrections 
are, at least at low redshifts, in reasonable agreement with 
those determined observationally by |Bell et al.| ( p003[ ). For 
example, in the K-band at z = 0.5 we find k and evolution 



(2003) 



Figure [12] shows that the space density goes up at the 
high luminosity end and down at the low luminosity end as 
the redshift decreases. These evolutionary changes have been 
quantified in Table [3] in terms of the Schechter parameters 
4>*/h 3 , M* — 51og 10 h and a. The model results indicate that 
the mass build-up of ellipticals is reflected in the if -band 
mostly through the large ~ 226% increase in <f>*/h 3 since z — 
1, and a small, 0.2±0.1 mag, brightening of M* — 51og 10 h. 

Parameterizing the LFs allows the luminosity density j 
to be calculated by integration, 



j = / 0(M)lO°- 4(M °- M) dM 



(3) 



where Mq is the absolute magnitude of the Sun, corre- 
sponding to 3.32 mag at K (Bell et al.|2003[ ). The luminosity 
density is important as it can provide a constraint on the 
mass density of stars in galaxies, given a mass to light ra- 
tio. The total if-band luminosity density at each redshift 
was calculated by integrating equation 3 over the interval 
—25 Mk— 51ogio/t ^ — 19 mag and the results are summa- 
rized in Table [3] These values indicate that the luminosity 
density, or equivalently, the mass density has increased by 
about a factor of 2.4 since z = 1. The predicted evolution 
since 2 = 2 is even more dramatic. These model results con- 
cur in general with recent observational results suggesting 
that the stellar mass contained within the red galaxy pop- 
ulation has at least doubled since z ~ 1 (|Bell et al.| 2004| 



Willmer et al.|2006| |Blanton|2006| fFaber et al.|2007||Brown| 



et 



al.|2007 l 



Our model results further indicate that the mass den- 
sity of luminous ellipticals grows at the expense of their 
lower luminosity counterparts causing the faint end slope 
to turn downwards with decreasing redshift. Such low lu- 
minosity ellipticals are presumably satellites of their more 
luminous counterparts, they are difficult to detect observa- 
tionally with current technology and contribute negligibly 
to the luminosity density evolution of ellliptical galaxies. 
At the other end of the luminosity range, recent observa- 



tions (Brown et al. 2007 Cool et al. 20081 indicate that 



the luminous Mb— 51ogio/i ^ -21 red galaxy population, 
loosely identified with ellipticals, were already in place at 
z = 1 whereas the corresponding model elliptical luminosity 
density, integrated between -24.5 Mk— 51ogio/i ^ -25.5, 
shows an increase by a factor of 4 since z = 1. Thus, the 
model predicts significant evolution of the most luminous 
elliptical galaxies since z = 1 that observers do not see. 
One possible reason for the discrepancy is that the colour 
selection criteria employed by observers may result in signifi- 
cant contamination of the observational datasets with galax- 
ies that are not actually ellipticals, causing elliptical galaxy 
evolution to be underestimated. This potential problem with 
the observations was noted previously by |Faber et al.| ( [2007[ ). 
It is crucial therefore to understand the connection between 
the luminous red galaxy population seen at z ~ 1 and the 
elliptical population at 2 = 0. Are the luminous red galaxies 
seen at 2 ~ 1 the progenitors of todays ellipticals or are they 
dusty objects destined to become what we recognise in the 
local universe as lenticulars, or early-type spiral galaxies? 
Of course, another possibility is that the model may be at 



The Origin of the Hubble Sequence in ACDM Cosmology 9 



Table 3. K-band Luminosity Function Fit Parameters 



Sample 


<j>*/h 3 


M»— 51ogi /i 


Q. 


j 




10 -4 galaxies Mpc~ 3 mag - 1 


mag 




1O 7 7iL Mpc" 3 


z = 


16.1 ± 0.5 


-23.93 ± 0.07 


— 0.48 ± 0.02 


10.4 


z = 1 


7.1 ± 0.3 


-23.77 ± 0.09 


-0.70 ± 0.02 


4.3 


z = 2 


4.9 ± 0.2 


-23.15 ± 0.08 


-0.82 ± 0.01 


1.9 



fault in overestimating the evolution of luminous ellipticals 
although, in defense of the model, it does succeed in pre- 
dicting the observed space density of ellipticals at z = (jj 
and the observed type averaged merger rate ( j |4.1[ ) and the 
observed star formation rate, as described in more detail in 
the next section. 



4.3 The Star Formation History as a Function of 
Morphological Type 

With regard to elliptical galaxies, the results presented 
in Figure [6] co nfirm the trend seen observationally (e.g., 
| Thomas et al.| ( |2005[ ) and references therein) that the star 
formation redshift for elliptical galaxies is correlated with 
their mass (luminosity) in the sense that the stars com- 
prising the more massive (luminous) ellipticals are old. The 
model predicts that half the stars comprising todays lumi- 
nous, —25^ Mr — 51og 10 h ^ — 24, ellipticals formed between 
redshifts of 2 and 4 with ages ~ 10 Gyr and were subse- 
quently assembled via mergers, that occurred since z ~3, 
into the objects we recognise in the local universe as lu- 
minous elliptical galaxies. In contrast, stars comprising the 
lower luminosity, — 24^Mk — 5 log 10 h, ellipticals are younger 
with ages spanning 4-10 Gyrs and that were created in star 
formation episodes the most recent of which were triggered 
by disk instabilities that occurred since z ~ 2. Thus, the 
seemingly anti-hierarchical behavior depicted in Figure[6]re- 
sults from the fact that there are different mechanisms gov- 
erning the star formation history of ellipticals: mergers for 
the luminous ones and disk instabilities for the less luminous 
ones. These results on formation redshift are insensitive to 
merger tree mass resolution (with galaxy mean stellar ages 
changing by only a few percent if we increase the merger tree 
resolution by a factor of 16) and to the details of the treat- 
ment of disk instabilities (adopting the minimalistic model 
of Appendix |A2] results in changes to mean stellar age of less 
than 10%). We note that, in our model, the amount of mass 
brought in to low luminosity ellipticals by minor mergers is 
relatively small due to the strong supernovae feedback which 
acts to suppress galaxy growth in lower mass halos. 

In Figure [13] we compare model predictions concerning 
the evolution of the specific star formation (SSF) rates with 
those measured recently by Perez-Gonzalez et al. ( 2008 1 for 
massive > 10 11 Mr. 



galaxies. The measured star formation 



7 Since the parameters of the model were adjusted to match the 
observed morphologically segregated luminosity functions this is 
not strictly a prediction. However, the unadjusted |Bower et al.| 
l |2006| ) model also yields a similar space density for elliptical 
galaxies as shown in Figure [T] 



rates are computed assuming a Chabrier initial mass func- 
tion and so correspond to the total star formation rate (i.e. 
including stars of all masses). The same is predicted by the 
model (in which we compute the total star formation rate 
from the known dynamical time and gas content of each 



model galaxy). The principal merits of the Perez-Gonzalez 
|et al.| p008| work are that the SSF rates are based on UV 
to Mid IR photometry and the results are segregated by vi- 
sual morphology. The figure shows that the observed SSF 
rates for 'disky galaxies' and 'spheroids' are bounded by the 
model predictions for 'disk' and 'bulge+disk' galaxies and 
both the model and observations show a well established 
trend that the SSFR is decreasing with decreasing redshift, 
with implications for the stellar mass accumulated in these 
types of systems ( |Dahlen et al.||2007| |Chen et al.|2009 l. 



4.4 The Secular Evolution of Bulges 

The model predicts that the bulges in spiral galaxies have 
grown through a combination of disk instabilities and minor 
mergers. Observationally, the colour of bulges are more sim- 
ilar to the disks in which they reside than to other bulges 



(Balcells & Peletier 1994 Peletier & Balcells 19961. This 



coupling between bulge and disk colour has been observed 
out to z~l (IDommguez-Palmero & Balcellsl [20081 [20091 



where younger, bluer high surface brightness bulges have 
been observed, suggesting an epoch of bulge formation. In 
concurrence with the observations, the model predicts that 
bulges in luminous A/k — 5 log 10 7l~ — 23 spiral galaxies have 
grown significantly through a combination of disk instabil- 
ities and minor mergers since z~2 and the model further 
predicts that bulges have been growing for a much longer 
time period, since z~4, in lower luminosity systems. 



5 CONCLUSIONS 

The Galform semi-analytic model of galaxy formation 
( Bower et al.|2006 1 has been employed to reproduce and ex- 
plain trends in galaxy morphology seen across a wide range 
of luminosities. Small adjustments to the model parameters 
permit excellent agreement with the observed morphologi- 
cally selected i^-band LF of 1 



the unmodified Bower et al. 



Devereux et al. (12009b , although 
( |2006[ ) model also performs very 



well. While predicting morphology from semi-analytic mod- 
els (and numerical simulations) remains challenging, a con- 
servative approach has been adopted by categorising model 
galaxies into just three broad classes: E ("pure spheroid"), 
SO-Sbc ("intermediate") and Sc-Scd ("pure disk") based on 
their K-band bulge-to-total luminosity ratio. Within this 
conservative categorisation, the results presented are con- 
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M K -51og(h) = -25.7 
B/T K =0.34 (50-Sbc) 



M K -5log(h) = -23.0 
B/T K =0.47 (SO-Sbc) 




0.5 



M/M n 



M/M n 



M K -51og(h) = -21.0 
B/T K =0.40 (SO-Sbc) 




M,(-51og(h) = -18.4 
B/T K =0.40 (SO-Sbc) 



1 0.5 1 

M/Mnn* M/M nnal 

tree diagrams for SO-Sbc galaxies shown for four different magnitudes as indicated in the panels. 
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Figure 11. The fraction of galaxies with mass >2.5 X 1O 1O M 
(corresponding to Mk — 51og 10 hi^. — 22.7) that have undergone 
a major merger (Af2/JWl>0.25) within the past 1 Gyr. The black 
line shows the predicted merger fraction averaged over all mor- 
phological types, while the blue, green and red lines show the 
predicted merger fraction for disk, disk+bulge and bulge morpho- 
logical classes respectively. Points show estimates of the merger 
fraction based on observations described by |Jogee et al.| ( |2009| ) 
and[L6pcz-Sanj uan et aTj j2009) (circles and star respectively). 



1 1 1 1 
_ z = 
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Figure 12. The predicted evolution of the elliptical galaxy lumi- 
nosity function with redshift. The luminosity functions have been 
k-corrected and passively evolved to z = 0. So, the differences 
reflect the mass build up, due primarily to mergers at the high 
luminosity end: -24.5 $C Mj^— 51ogio/i $J -25.5 



verged with respect to the resolution of the N-body-derived 



merger trees (see Appendix All. Our conclusions follow 



• The Sc-Scd disk galaxies, with ~ 10% or less of their 
K-band light coming from the bulge, can form in abundance 
even in a hierarchical cold dark matter cosmology in which 
merging is, in general, commonplace. The median merger 
history of the halos that these galaxies inhabit is not signif- 



■s 




Figure 13. The evolution of the specific star formation rate with 
redshift for massive 10 11 Mq galaxies. The blue and red crosses 
identify the disky galaxies and spheroids of Perez-Gonzalez et al. 
(2008). The blue and green lines identify the model disk and 



disk+bulge galaxies respectively. The model SSF rates for the 
bulge dominated elliptical galaxies are essentially zero. 



icantly different from that of elliptical galaxies which do ex- 
perience significant merging. However, the variation around 
that median history is much less for the Sc-Scd galaxies 
than for ellipticals. The Sc-Scd class typically experience a 
very simple accretion history, deviating little from the me- 
dian and exhibiting relatively few large merging events. The 
formation history of these systems is therefore determined 
by the rate of growth of their dark matter halo and the 
corresponding rate of gas supply, cooling and infall. Galax- 
ies in this class have typically been forming stars in their 
disks over large fractions of cosmic history and continue to 
do so at the present day. Steinmetz & Navarro ( 2002| and, 
more recently, |Dutton| ( |2009| | have suggested that produc- 
ing bulge-less disk galaxies in sufficient numbers could be 
challenging in a CDM universe. In fact, because there is a 
sufficiently large number of dark matter halos with relatively 
quiet merging histories, bulge-less disk galaxies can prevail, 
in the numbers observed. Additionally, strong feedback pre- 
vents the formation of relatively massive satellites around 
these galaxies which are predominantly low mass. 

• Luminous elliptical galaxies form through multiple 
mergers of smaller progenitors. For the most massive ellipti- 
cals these mergers tend to be "dry" and so contribute no new 
star formation. Interestingly, the Galform model does not 
produce the dwarf elliptical sequence which should emerge 
at Mk — 51og 10 /i ~ —21 mag. Thus, the model prediction 
concerning the shape of the LF for low luminosity ellipticals 
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(with Mk — 51og 10 h > —21 mag) and the model prediction 
that they form through disk instability events is uncertain. 

• The intermediate morphological types, SO-Sbc, are 
characterized by an early period of merging and disk insta- 
bility events which form the bulge, followed by an extended 
period of steady, uninterrupted disk growth. 

• The model results concur with several observational re- 
sults including the type averaged galaxy merger rate, mea- 
surements of the specific star formation rate, the secular evo- 
lution of spiral bulges, the seemingly anti-heirarchical evolu- 
tion of elliptical galaxies and the factor of ~ 2 growth in the 
luminosity density of elliptical galaxies since z ~ 1. With re- 
gard to the last point however, the model predicts significant 
evolution at the high luminosity end of the A'-band elliptical 
galaxy LF that has yet to be confirmed observationally. 
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Figure Al. A study of resolution effects on the morphologically 
segregated luminosity function of Fig. [I] The thickest lines show 
the results obtained utilizing Millennium Simulation merger trees 
(i.e. as used throughout our work). The remaining three sets of 
lines were produced using analytically generated merger trees (us- 
ing the algorithm of |Parkinson et aT] ( |2008"| ). The thickest of these 
lines correspond to analytic trees with a mass resolution matched 
to that of the Millennium Simulation, while successively thinner 
lines have 4 and 16 times better mass resolution. 



APPENDIX A: ROBUSTNESS OF MODEL 
MORPHOLOGIES 

In this appendix we explore to effects which could influence 
the morphological properties of model galaxies, and demon- 
strate that, in fact, our results are robust to these effects. 



Al Resolution Study 

Since morphological evolution in our model is driven at 
least partly by merging activity, it is important to assess 
whether the resolution of the N-body merger trees that we 
utilize causes us to miss some lower mass merging activity 
and, therefore, to incorrectly estimate morphologies of some 
galaxies. While we currently do not have suitable N-body 
merger trees with higher resolution available to us for a reso- 
lution study, we can instead utilize merger trees constructed 
using a modified extended Press-Schechter approach (see 



Parkinson et al. 2008 1 which is designed to produce trees 
which are statstically equivalent to those found in N-body 
simulations. The advantage of this approach is that it allows 
us to increase the resolution to almost arbitrarily high lev- 
els. While the results will not precisely match those obtained 
using N-body trees they are sufficiently close to allow us to 
explore the effects of changing resolution. Figure |A1| below 
shows the results of a resolution study carried out in this 
way on the morphologically segregated luminosity function. 

The thickest lines correspond to the case of N-body 
merger trees drawn from the Millennium Simulation (i.e. 
used in this paper). The remaining lines were all generated 



using the exact same model parameters but with merger 
trees constructed using the |Parkinson et ah| ( |2008| l algo- 
rithm. Of these lines, the thickest has a resolution matched 
to that of the Millennium Simulation. As can be seen 



(Fig. Al I, it does not produce luminosity functions in precise 
agreement with those obtained using Millennium Simulation 
merger trees (since the analytic merger trees are not exactly 
statistically equivalent to the N-body ones), but they are 
very close. The thinner lines represent results obtained when 
the mass resolution in the analytic merger trees is increased 
(i.e. shifted to smaller masses) by factors of 4 and 16. It 
can clearly be seen that there is no significant change in the 
predicted luminosity functions. As such, we conclude that 
galaxy morphologies are well converged in our model when 
using trees with the Millennium mass resolution. (This is 
a consequence of the effects of supernovae feedback, which 
makes galaxy formation in lower mass halos highly ineffi- 
cient and, therefore, to have little effect on more massive 
galaxies.) 



A2 Alternative Treatment of Disk Instabilties 

The instability of galaxy disks is a difficult process to treat 
analytically. Our standard approach assumes that any un- 
stable disk will quickly develop a strong bar and destroy 
itself entirely, turning the galaxy into a pure spheroid. Re- 



cent, high resolution hydrodynamical simulations (Agcrtz 
|et al.|2009[ ) show that unstable disks experience local insta- 
bilities leading to fragmentation and the incomplete destruc- 
tion of the disk. To at least partially assess the importance 
of the assumptions we make regarding disk instabilities we 
have implemented an alternative disk instability treatment 



16 Andrew J. Benson and Nick Devereux 




-20 -22 -24 

M K- 5l °eio h 



Figure A2. The effects of an alternative disk instability model on 
the morphologically segregated luminosity function. The thickest 
lines show results from our standard treatment of disk instabili- 
ties applied to merger trees extracted from the Millennium Simu- 
lation. Intermediate lines indicate the same standard treatment of 
instabilities but applied to analytically derived merger trees (us- 
ing the algorithm of |Parkinson et ah] ( |2008^ ). Finally, the thinnest 
lines show the results of adopting an alternative treatment of disk 
instabilities (in which only just enough mass is transferred from 
an unstable disk to cause it to restabilize) applied to analytically 
derived merger trees. 



in our model in which only just enough mass is transferred 
from the disk to the spheroid component to restabilize the 
disk — an approach which has been adopted by other semi- 
analytic models ( Hatton et al.|2003 Croton et al. 2006). Fig- 
ure |A2| shows the effect of switching to this alternative treat- 
ment on the morphologically segregated luminosity function. 
The thickest lines show results utilizing Millennium Simula- 
tion trees with our standard treatment of disk instabilities. 
Medium thickness lines show the same treatment of instabil- 
ities but with analytic merger trees with resolution matched 
to that of the Millennium Simulation. Finally, the thinnest 
lines show results of using the alternative treatment of in- 
stability events with analytic merger trees. 

The net result of switching to the alternative treatment 
of instabilities is to slightly increase the number of bulgeless 
galaxies at all luminosities, with a corresponding decrease 
in the numbers of intermediate and pure spheroid galaxies. 
The changes do not alter the qualitative trends of morpho- 
logical mix with luminosity. However, in order to reproduce 
quantitatively the observed morphological mix of galaxies 
at Mk — 51og/i = —23 would require adjusting the range of 
B/T ratios used to segregate galaxies into the three broad 
morphological classes (Table [T| . The sense of the required 
adjustment for the E ("pure spheroid") class would be to 
smaller B/T ratios, which is inconsistent with their observed 
morphology and may be a sign that our standard treatment 
of disk instabilities is more realistic. 



